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Abstract 

Spread spectrum, widely employed in modern digital wireless terrestrial radio systems, chooses 
a signal with a noise-like character and much higher bandwidth than necessary. This paper ad- 
vocates spread spectrum modulation for interstellar communication, motivated by robust immu- 
nity to radio-frequency interference (RFI) of technological origin in the vicinity of the receiver 
while preserving full detection sensitivity in the presence of natural sources of noise. Receiver 
design for noise immunity alone provides no basis for choosing a signal with any specific charac- 
ter, therefore failing to reduce ambiguity. By adding RFI to noise immunity as a design objective, 
the conjunction of choice of signal (by the transmitter) together with optimum detection for noise 
immunity (in the receiver) leads through simple probabilistic argument to the conclusion that the 
signal should possess the statistical properties of a burst of white noise, and also have a large 
time-bandwidth product. Thus spread spectrum also provides an implicit coordination between 
transmitter and receiver by reducing the ambiguity as to the signal character. This strategy re- 
quires the receiver to guess the specific noise-like signal, and it is contended that this is feasible 
if an appropriate pseudorandom signal is generated algorithmically. For example, conceptually 
simple algorithms like the binary expansion of common irrational numbers like tt are shown to be 
suitable. Due to its deliberately wider bandwidth, spread spectrum is more susceptible to disper- 
sion and distortion in propagation through the interstellar medium, desirably reducing ambiguity 
in parameters like bandwidth and carrier frequency. This suggests a promising new direction in 
interstellar communication using spread spectrum modulation techniques. 
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1. INTRODUCTION 

This paper addresses the design of an end-to-end digital interstellar communication system that 
can be used to exchange information among intelligent civilizations. As advocated in [IJ, we 
employ communication's engineering principles to achieve an implicit form of coordination. This 
strategy leverages the physical impairments that are unavoidable in an end-to-end interstellar dig- 
ital communication system, and presumably known to both transmitter and receiver, and bases 
design decisions on mitigation of those impairments. Where those design decisions are based on 
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mathematical optimization, the transmitter and receiver designers will necessarily arrive at com- 
patible designs. Where optimization is not possible, the receiver must apply judgement and/or 
systematically search over the possibilities, minimizing ambiguity by keeping the design as simple 
as possible within the constraints of acceptable cost and performance. 

There are two distinct phases to establishing information transfer lH: (a) receiver discovery of the 
existence of a signal followed by (b) extraction of information from that signal. This is rendered 
more challenging by the impossibility of advance coordination between transmitter and receiver, 
and discovery is the more challenging problem and also represents the dominant cost and effort 
for a receiving civilization. Although this paper addresses both discovery and communication 
from an end-to-end perspective, past efforts have addressed transmission and reception as pri- 
marily separate issues, and in the case of the receiver (with the exception of [2 , 3|) have focused on 
discovery of signals designed to attract attention but not bear information. Some entryways into 
the extensive literature on this topic can be found in ||4l|5l|6j[Zl|8l|9l. Longstanding programs try 
to discover signals of technological origin [10 1 (called "SETI" fi^). Other efforts have transmitted 
signals with an embedded message (called "METI" [llj), including the "Arecibo message" [12J in 
1974 and more recently "Cosmic Call", "Teen Age Message" and "A Message from Earth" fSl. The 
transmitter side has also been studied from a cost perspective Il3tll4l . 

The end-to-end design of interstellar communication systems is relevant to both SETI and METI. 
Our focus is on radio frequencies, digital modulation (although analog modulation has also been 
experimented with in METI), and communication of an information-bearing message [9] . Many 
SETI projects have followed the Cyclops Report |15| more than four decades ago, which recom- 
mended searching for extremely narrow bandwidth signals. Several factors have changed since 
then. First, the explosion in terrestrial wireless multi-user communications systems has improved 
our understanding of radio communications in an RFI-domrnated environment IIT6l[T7l[T8l . This 
has strongly shifted the design of radio systems toward spread spectrum techniques, in which 
the signal bandwidth is much greater than dictated by the supported information rates |19 , 20| . 
Second, there has been great progress in understanding the radio propagation properties of the 
interstellar medium (ISM) based in large part on pulsar astronomy [21J. Third, dramatic advances 
in electronics technology open up new possibilities in our own search capabilities. 

Here we consider the specific design criterion of maximizing the probability of discovering an 
information-bearing signal and the correct detection of its embedded information content. This 
criterion is appropriate for a transmitter constrained in transmit power and seeking to maximize 
the distance over which a message signal can be discovered and its information content extracted. 
For a given transmit power this criterion also maximizes the number of stellar systems wherein 
discovery is feasible, and in this sense the probability that the signal will be discovered by some- 
one somewhere. Even if the details of the transmit signal are completely known to the receiver, 
there are four primary impairments beyond the control of transmitter and receiver designers that 
impinge on discovery and detection: natural noise sources, radio-frequency interference (RFI) in 
the receiver vicinity, radio propagation effects [22] through the interstellar medium (ISM), and 
Doppler due to relative accelerations. This paper focuses on the first two, and subsequent papers 
(starting with Il23l ) will address the remaining. 

RFI is a challenge for SETI because its technological origin mimics the interstellar signal being 
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sought. An ongoing trend on earth and its immediate environs is ever-increasing local RFI. It is 
possible that a transmitter designer (which is likely to be more technologically mature than us) has 
experienced a similar growth in RFI and also has developed an understanding of RFI mitigation 
in the context of its own multi-user wireless communication systems. If so, the transmitter de- 
signer may view RFI mitigation in the receiver environment as a significant design objective. We 
demonstrate in this paper that this circumstance would be very fortunate, for two reasons. First, 
signal choice by the transmitter taking RFI into account can strongly mitigate the adverse effects 
of RFI on our own SETI observations, directly addressing our growing challenge. Second, consid- 
erations of signal acquisition in the presence of RFI allow us to infer very specific characteristics 
of an information-bearing signal, considerably beyond what design for noise (as pursued in [15 J 
and elsewhere) can infer. This establishes principles allowing us to infer what type of signal may 
have been chosen by the transmitter. In contrast, even recent attempts at applying communication 
theory to SETI have emphasized noise sensitivity [24T l25ll26l l2ll27ll28l, and are as a result unable 
to infer any concrete characteristics of the signal even as many advocate the use of greater band- 
width. The propagation characteristics of the ISM introduce further constraints on the bandwidth 
and carrier frequency of the transmit signal, which desirably limits the ambiguity in these param- 
eters for both transmitter and receiver 1.23 J . Others have emphasized such astrophysical effects in 
trying to infer signal characteristics ||29| . 

The focus of this paper is on the choice of a transmitted signal, which directly parallels the re- 
ceiver's challenge of anticipating what type of signal to expect. In this we take the perspective 
of a transmitter designer, because in the absence of explicit coordination it is the transmitter, and 
the transmitter alone, that chooses the signal. This is significant because the transmitter designer 
possesses far less information about the receiver's environment than the receiver designer, due 
to both distance (tens to hundreds of light-years) and speed-of-light delay (tens to hundreds of 
years). While the receiver design can and should take into account all relevant characteristics of its 
local environs and available resources and technology, in terms of the narrower issue of what type 
of signal to expect the receiver designer must rely exclusively on the perspective of the transmitter 
designer. 

The primary result of this paper is as follows. The transmitter designer should assume that the 
receiver employes the optimum detection algorithm for white Gaussian noise, but assume noth- 
ing specific about the RFI environment of the receiver since nothing is presumably known. The 
transmitter can, however, explicitly design a transmit signal that minimizes the effect of RFI on 
the receiver's discovery and detection probabilities in a robust way; that is, in a way that provides 
a constant immunity regardless of the nature of the RFI. It is shown that the resulting immunity 
increases with the product of time duration and bandwidth, and that the signal should resem- 
ble statistically a burst of white noise. Intuitively this is advantageous because RFI resembles 
such a signal with a likelihood that decreases exponentially with time-bandwidth product. Both a 
transmitter and receiver designer using this optimization criterion and employing the tools of el- 
ementary probability theory will arrive at this same conclusion. Although the context is different, 
variations on this principle inform the design of many modern widely deployed terrestrial digital 
wireless communication systems, so this has been extensively tested in practice and is likely to 
have a prominent place in the technology portfolio of an extraterrestrial civilization as well. 
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2. DESIGN APPROACH 



This section qualitatively describes an end-to-end digital communication system including coun- 
teracting noise and RFl. Section |3] brings to bear a relevant mathematical principle drawn from 
probability theory. 



2.1. Modulation 



An intermediate step of modulation I Appendix A| accommodates the analog nature of a ra- 



dio channel. A common case is pulse-amplitude modulation (PAM), in which a sequence of 
information-bearing symbols {Aj.} , each drawn from a finite constellation (finite set of discrete 
points), is used to adjust the amplitude of a sequence of pulses with shape h{t), 

y{t) = ^h-E^kh{t-kT). (1) 

k 

It is assumed that h{t) is a unit-energy pulse waveform, 1/T is the symbol rate, and f;, is the 
energy per symbol. If the receiver can distinguish which Aj. was transmitted, it can extract the 
information being conveyed. The average power is Vs = £s/ T, and if it is desired to reduce Vs (at 
the expense of lower information rate) this can be done by increasing T. The radio propagation 



must be passband with some carrier frequency fc, and as described in Appendix A this allows the 
baseband-equivalent signal y{t) in to be complex-valuedj^ This paper focuses on the choice of 
h{t) and its implications to the extraction of {Aj-} from i/{t) in the face of noise and RFl. 

2.2. Detection in noise 

Thermal noise is ubiquitous in nature; at radio frequencies, it can be modeled at baseband as an 
additive white Gaussian noise N{t). Here we focus on the detection of representative one symbol 
Aq. lfh{t) is known to (or can be guessed by) the receiver, the optimum receiver processing 1,30.1 is 
a matched filter (MF), 

Z = J {yit) + N{t)) h*{t)dt = ^sAo + J N{t)h*{t)dt. (2) 
For ||2} to be valid the time translates of h{t) must also be orthogonal, 
J h{t - kT) h* {t) dt = or J2 ^ {f + ' 



(3) 



The decision variable Z in ||2} can be used for two purposes. During commimication, Z can be 
applied to a complex-valued slicer that takes account of the symbol constellation to make a de- 
cision on which Aq was transmitted. During discovery the symbol constellation is not known, 
but as long as h{t) is known or can be guessed |Z| (an estimate of \Ao\) can be applied to a 



^ At passband A^. is interpretated as the carrier amplitude and phase for symbol period k. 
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real-valued threshold to decide whether or not a signal is present. Thus, the MF can be used for 
discovery based on single-symbol energy, the approach used here. If more assumptions are made 
about the symbol constellation and any redundancy in successive symbols, then the approaches 
here can be modified to base discovery on multiple symbols, thereby capturing more signal energy 
(see Is) for an example). 

The statistics of the noise term in Z in ||2} do not depend on h{t), including its time duration or 
bandwidth]^ As a result, design issues related to noise offer no clue as to what h{t) should be 
chosen by the transmitter. This is paradoxical since the MF also assumes that the receiver knows 
or can guess h{t). 

2.3. Radio-frequency interfrrence 

Radio astronomy is increasingly impacted by RFI with increasing receiver sensitivity and grow- 
ing sources of RFI [ST/S?!, and many techniques have been successfully applied to RFI mitigation 
133] . This growing problem has even stimulated extreme measures like planning for observations 
on the far side of the moon ||34ll35| . Many of the same instruments, challenges, and techniques 
can be shared with interstellar communication, including those listed in Table [l] Interstellar com- 
munication has the special challenge that RFI and the interstellar communication signal share 
a technological origin, and may even be communication signals using similar modulation tech- 
niques. In the discovery phase, RFI creates two very different anomalies: False alarms due to 
RFI being mistaken for signal, as well as misses due to RFI overpowering a signal of interest. A 
special opportunity applicable only to communication (technique V. in Table [TJ is to engineer the 
signal structure to make it distinguishable from RFI. This approach is familiar in terrestrial and 
space-based communication systems, since they also have to deal with RFI and the challenge of 
distinguishing one radio-based communication from another. Gains arising from this approach are 
additive to gains from other mitigation measures, and unlike many techniques used in astronomy 
they can be effective against RFI with time, spectral, and spatial overlap with the communication 
signal, reducing both false alarms and misses. A disadvantage is the required cooperation of the 
transmitter. 

This paper is dedicated to investigating this signal structure approach. To address the effect of RFI 
on communication, observe that the matched filter of l|2} does not take account of RFI. If further 
knowledge of the characteristic of RFI were incorporated into the detection, then sensitivity could 
be improved. However, since the transmitter designer presumably has no specific knowledge of 
RFI local to the receiver, that designer should choose h{t) assuming only that the noise-optimal MF 
is used for detection. Even subject to this limitation, the transmitter designer can assist the receiver 
through the choice of h{t), and that observation provides specific and credible information about 
the nature ofh{t) that a receiver should seek. 

Suppose the baseband-equivalent RFI at the receiver is r(f). The transmitter designer should take 
the stance that r ( ■ ) is an unknown waveform, and even be unwilling to assign a statistical distri- 
bution (which assumes that some r( )'s occur more frequently than others). The transmitter can 



^ The statement in (6) that "spreading the bandwidth makes detection more difficult" does not apply to a MF. 
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Table 1: Some techniques for dealing with RFl borrowed from radio astronomy and their efficacy for inter- 
stellar commimication. It is assumed the receiving antenna is highly directional and that RFI originates from 
a point source. The terminology is DOA = "direction of arrival" and SOI = "signal of interest". 



Technique 



Interstellar communication 



I. Spatial 
rejection 



II. Temporal 
rejection 



III. Frequency 
rejection 



IV. Spatial 
filtering 



V. Signal 

character 



Situate receiving antennas in quiet loca- 
tions. Reduce spurious antenna side- 
lobes by antenna design. Space-based RFI 
sources not in the SOI's DOA are rejected. 

If an RFI source is transient, observe dur- 
ing periods it is absent. Most space-based 
sources of RFI (even geostationary satel- 
lites) are in motion relative to a SOI's DOA. 
RFI may exhibit distinguishing periodici- 
ties |36|. 



Effective, although an RFI source in the 
DOA is difficult to distinguish from a SOI 
since they are both of technological origin. 

The SOI will likely be transient (due to stel- 
lar scanning, random scintillation, etc.), but 
its periodicity can be chosen to be non- 
typical of near-space RFI. 



If an RFI source is frequency-specific, ob- 
serve at other frequencies. RFI may have 
a recognizable narrowband component sig- 
nature. Some frequency bands are reserved 
for passive use (no transmissions). 



A discovery search cycles through fre- 
quency, and only RFI overlapping the 
SOI bandwidth can cause false alarms or 
misses. Passive frequency bands offer less 
protection since a remote transmitter is not 
cognizant of their location. 



For multiple-element antennas, take ad- 
vantage of high or low correlation of RFI 
across elements and different DOA's of RFI 
and SOI L37J . Employ array signal process- 
ing, such as RFI tracking, adaptive interfer- 
ence cancellation, etc. 



Effective, with the same caution as in I. 



SOI's may have a special structure (like 
the periodicity of pulsars |38|) that distin- 
guishes them from RFI or makes interme- 
diate ISM impairments evident Ii39i . 



With the cooperation of the transmitter, en- 
gineer the SOI's structure to make it more 
readily distinguishable from RFI. 



anticipate the magnitude of the contribution of RFI to Z in ||2} as 



r{t)h*(t)dt 



R{f)H*{f)df 



<£r{T,B), 



(4) 



where it is assumed that h{t) is simultaneously bandlimited to / E [0,B] Hz (so the passband 
representation is bandlimited to |/| £ [fcfc + B]) and time-limited to t E [0, T] secondsPlThis RFI 
term corrupting the decision can be small or large in magnitude depending on how r ( ■ ) correlates 



^ This arbitrary choice of fc (so that the baseband signal has a one-sided spectrum) emphasizes that h{t) is complex- 
valued. The ideal of simultaneous time and bandwidth limitation can be achieved with increasing accuracy |40| as -K = 
BT oo. 



6 



(or doesn't correlate) with signal h{t), but its largest value is £r{T, B), the energy of r{t) confined 
tote [0,T] and / e [0,B]. 

With respect to RFI, the receiver gain G > 1 is defined as the ratio of the input RFI energy to the RFI 
energy appearing at the decision threshold. The larger G the better the receiver has suppressed RFI 
before decision making. From Q there are two obvious components in G = FG x PG. The filtering 
gain FG = £,.{00,00)/ £,.{T,B) is the reduction in energy due to confinement to i e [0, T] and 
/ E [0, B] (an example of techniques II. and III. in Table[l|. The processing gain PG (an example of IV. 
and V.) is the reduction in energy below £r{T, B) due to a less-than-maximum correlation between 
r(f) (after confinement in bandwidth and time) andh{t). Clearly the worst case is r(t) — y/£^h{t), 
for which G = FG = PG = 1 . Not surprisingly, a receiver cannot distinguish h{t) from an identical 
RFI (within a constant amplitude and phase shift). 

3. ISOTROPIC PRINCIPLE AND ITS CONSEQUENCES 

An appropriate choice of h{t) can insure that PG ^ 1 consistently in a statistical sense. One key 
is to choose K = BT ^ 1. This is called spread spectrum because B is larger than the minimum 
required, which is B > 1/T from ||3| (this is called the Nyquist criterion |30|). Thus large PG 
results in part from choosing a bandwidth B much larger than the minimum required for a given 
symbol rate T. This will follow from the isotropic principle (IP), which explains both the optimality 
of the MF in Q and the relationship between bandwidth and PG. 

3.1. Signal degrees of freedom 

Assume that h{t) equals a weighted linear combination of K = BT orthonormal functions, imply- 
ing that h{t) actually falls in a K dimensional space. Common examples include 

1 ^^^/k + l/2\ ft-k-l/2\ 
M0 = ^ E^^ic^''"'''" (6) 

k=0 

The first form in ||5} is the sampling theorem, which uses time translates of a appropriate imit-area 
interpolation function <^() to generate an that is band limited to / e [0,B] and approximately 
time limited to f G [0, T]. The second form in ||6} is the Fourier series, which generates an h{t) that 
is time-limited to i G [0, T] (with periodic extension) and approximately band limited to / G [0, B] . 
Both expansions illustrate that B and T can be chosen independently (as long as BT > 1), and this 
will not affect the detection noise immunity as long as £s is fixed [^Increasing T does reduce both 
Vs and information rate. 



* The vague notion of waveforms that are both band limited and time limited can be made precise by choosing a set of 
orthonormal functions as these examples illustrate and then defining a linear span of these functions |40|. The assertion 
in 1 15 1 that increasing T in the interest of reducing signal power "narrows the signal spectrum" and as a consequence "we 
would expect interstellar contact signals to be highly monochromatic" is unjustified because B can be chosen independently 
of T. 
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Table 2: Three flavors of complete randomness (CR) for a random V 6 C . 



White CR All 2K real and imaginary parts are zero-mean, mutually uncorrelated, and 

have the same variance /2. 

Independent CR All 2K real and imaginary parts are zero-mean and statistically independent. 

Since independent random variables are uncorrelated, independent CR implies 
white CR. 

Gaussian CR White CR, and in addition all 2K real and imaginary parts have a joint Gaussian 

distribution. Since uncorrelated Gaussian variables are independent, Gaussian 
CR implies independent CR. 



An orthonormal representation converts the problem of designing h{t) to the choice of a h G 
(Euclidean space of X-dimensional complex-valued column vectors). With respect to any basis, 
the reception after band limiting to / e [fc, fc + B] and translation to baseband / E [0, B] can be 
represented by a random Y G as 

Y= y^e''^-h + Ni + v/^-r, (7) 

where £)j is an unknown signal energy, is an unknown carrier phase, h is a known unit-energy 
pulse, Nj is a random noise vector, £r is an unknown RFI energy, and r is an unknown unit-energy 
RFI vector. The energy of any x G is ||x|| ^ = x^x, where x^ denotes the transpose and conjugate 
of X. 

Any h points in a single direction in C*', and K specifies the number of complex numbers required 
to specify that direction, or equivalently the dimensionality of the space containing h. Thus K, 
which we call the signal degrees of freedom (DOF), quantifies the ambiguity in the signal direction. 
Even h = [1, 0, . . . , 0]^ has a DOF of K; it just happens that K — 1 of those degrees have been 
chosen to be zero. 

3.2. Completely random vectors and the isotropic principle 

Completely random (CR) vectors play a fundamental role in both noise analysis and signal design. A 
random V G is fully specified by 2_K real-valued random variables {5R{Vjt}/ ^'{Vj;}/ 1 ^ ^ < 
with some joint probability distribution. A Gaussian CR V is zero-mean (£ [V] = 0), and all of its 
2K constituent random variables are jointly Gaussian, identically distributed with variance 
and statistically independent. These conditions can be relaxed as shown in Table |2]by eliminating 
the Gaussian and/or independence conditions. 

CR vectors obey an isotropic principle (IP) | |4T] that comes in two flavors listed in Table |3] Roughly 
speaking, the IP says that statistics of a white CR vector are the same regardless of what direction 
in C^ we look and regardless of what orthogonal basis we choose. Thus, the energy (variance) of 
a white CR vector is spread uniformly across coordinates, and this remains true after any unitary 
coordinate transformation (such as the Fourier transform). 
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Table 3: Two flavors of the isotropic principle (IP) for a rando m V 6 C^, w here u is an arbitrary unit vector 



u|| = 1) and U is an arbitrary unitary matrix (UU = I). See Appendix B 



for derivations and elaboration. 



Energy 

isotropic 

principle 

(V is white CR) 



The component in any direction u has the same energy (E|u^V| ^ = cr^). 
V after coordinate transformation (UV) is also white CR. 



Distribution 

isotropic 

principle 

(V is Gaussian CR) 



u'''V is Gaussian, and thus 2 ■ |u'''V| ^/c7^ has a x'^{%) distribution. 
UV is Gaussian CR. 

If V is merely independent CR, then under mild conditions the distri- 
bution of u^V approaches Gaussian as — >^ oo. This implies that the 
distribution of 2 ■ \ I a'^ approaches a ;)(;^(2) distribution in the limit 

as — >■ 00. 



3.3. Noise 

At radio frequencies it is usually assumed that the passband noise from natural sources (e.g. the 
background radiation of the universe and thermal noise in the recei ver) is white a nd Gaussian, 
in which case Ni is Gaussian CR with variance a'^ per component (Appendix Cl. Our goal in 
detector design is to process Y to yield a real-valued decision variable Z which is applied to a 
threshold to distinguish > (signal plus noise and RFI) from f;, = (noise and RFI alone). 
Signal is declared present when Z > A, define a normalized threshold t = \/ a, and let PpA be the 
false alarm probability and be the detection probability. 



For the moment assume the absence of RFI (r = 0) and contrast two cases. 



h^Y = 



h+ Ni 



or 



Yll = 



■Ni 



(8) 



where the MF Zm calculates the magnitude of the component of reception Y in the direction of 
unit vector h, and the energy detector (ED) Z^ calculates the square root of the total energy in Y. 
The distribution Zm and Zg are determined in Appendix B The detector sensitivity is defined as 



the 8h required to maintain fixed PpA ai^d Pq as K and h change. 

For a given K, the sensitivity of neither the MF or the ED depends on the choice of h. It is in 
the dependence f;, on K that the MF and ED differ markedly. For K = 1 they are the same, and 
for X > 1 the sensitivity of the MF does not change. This follows immediately from the IP, since 
the distribution of Zm does not depend on either h or K. On the other hand, for the ED S^^ must 
increase by about 10 dB for each 10 x increase in K. This behavior is explained by Figure [ij where 
T must be increased by that amount to maintain fixed PpA due to the increasing total noise in Zg. 
This increase in t has to be compensated by an increase in by that amount to maintain fixed 
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loglO^FA 




Tin dB 



Figure 1: The energy detector log-log plot of PpA against the normalized threshold t for K = 10" where 
n 6 [0,4]. Larger K have higher PpA for the same t, or require larger t for the same PpA- Note that the 
matched filter corresponds to K = 1, n = 0, and has the most favorable tradeoff of t vs PpA- 



A major advantage of the ED is that unlike the MF it does not require knowledge of h, but a 
substantial penalty is paid in sensitivity for K > 1. That is one manifestation of a more general 
principle of signal detection: The more specific the detector's knowledge of the signal, the better 
the sensitivity of a detector that makes use of that knowledge. Cyclops IITSlI prescribed an ED to 
avoid knowing h, but also prescribed X « 1 in order to preserve detector sensitivity. Intuitively 
the penalty for increasing K = BT in the ED arises because more total noise energy is captured 
in Ze either because B is large (noise energy equals noise density times bandwidth) or T is large 
(noise energy equals noise power times integration time). K ~ 1 can be accomplished with a 
carrier-like signal (small B ~ 1/T) or a pulse-like signal (small T ~ 1/B). In both cases the one 
DOF specifies the phase and amplitude, and thus this strategy chooses an exact signal waveform 
h{t) (just as we propose for spread spectrum in Section|4j. Small B was preferred in 1 15 1 to reduce 
the transmit power for a given £}^. 



3.4. Radio-frequency interference 

Is there any motivation to increase the signal DOF K? First, the information-bearing nature of 
the PAM signal of Q is not a motivation, since the only requirement is that ^ be satisfied and 
this only requires that K = B T > 1 . Decreasing T is beneficial because the higher symbol rate 
implies (all else equal) a higher information rate, and this in turn requires larger B but not larger 
K. The primary penalty for increasing the information rate is higher transmit signal power, since 



that power is f / T, and ifh{t) is known and a MF is used in the receiver we saw in Section 3.3 that 
£h does not depend on B, T, or K for fixed PpA and Pp. There is, however, one strong motivation 
for making K large, and that is RFI as we will now see. 



Recall from Section 2.3 that the receiver's immunity to RFI is quantified by the gain G = FG x PG. 
If, as is likely, the transmitter designer has no specific knowledge of the RFI environment at the 
receiver, he has no opportunity to manipulate FG through the choice of T and B since the temporal 
and bandwidth characteristics of the RFI are unknown. However, that designer does have an 
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opportunity to manipulate PG by choosing K to be large. As long as the receiver has knowledge of 
h{t) and uses that knowledge to construct a MF detector, then the transmitter designer can assume 
that such a large K has no adverse effect on the detection noise sensitivity. He can therefore focus 
on choosing h{t) to achieve robust PG, meaning the greatest PG that can be consistently achieved 
without needing to incorporate any specific assumptions about the RFI. 

From the transmitter designer perspective, assume the RFI after the benefit of FG is represented 
by a deterministic (but unknown) vector r G in the same basis as the signal vector h, where 
r is a imit vector that represents the direction of the RFI in and is the energy of the RFI. If 
the receiver is presumed to use a MF, then the transmitter designer knows that the component of 
RFI at the decision threshold is Zm = £r h^r , and this term biases the decision in favor of > 
and hence increases PpA- Iri choosing h the transmitter thus has an interest in making | h^r | as 
consistently small as possible. The factor by which RFI energy is reduced by the MF is 

PG = II ^'•'^11^ = I (9) 

PG depends on the geometric relationship of two unit vectors, h and r: PG = 1 when h is colinear 
with r, PG = 00 when it is orthogonal, and otherwise 1 < PG < oo. 

How do we achieve consistently large PG? Two perspectives arrive at similar conclusions. First, 
if h is in some sense spread uniformly throughout K dimensional space, only a fraction 1/K of its 
unit energy will fall in the direction of any r, independent of r. In that case, PG w K, and there is a 
clear benefit to making K as large as possible. Second, consider a basis chosen to include h as the 
first basis vector Then two extreme scenarios are the unit energy of r spread evenly over all basis 
vectors, in which case only a fraction 1/K of the unit energy is co-linear with the one coordinate 
h, or the unit energy is concentrated in one basis vector, in which case the chance that this vector 
happens to be the signal direction is only 1/K ||T9| . In either case, PG w X in a statistical sense, 
and we conclude again that large K is advantageous. Nevertheless, PG ~ 1 can always occur 
inadvertently when h w r, even if that outcome is deemed to be unlikely. 

Lacking control over or knowledge of r, a communication engineer focuses instead on h, which 
the transmitter designer chooses. Specifically, the previous argument suggests that he should seek 
a signal which spreads its energy uniformly throughout the K DOF's. There is no direct justifica- 
tion for choosing any specific h, either for noise sensitivity or RFI immunity. Thus, the designer 
can do no better than choose h as a random vector H drawn from some appropriate random en- 
semble. The transmitter will choose H randomly from this ensemble, and a receiver knowing 
H implements a MF using H. The obvious question — how does the receiver know H? — is ad- 



dressed in Section 4.2 With this approach, the designer can answer two useful questions: First, 
PG = 1/ |H'''r| 2 becomes a random variable, and the distribution of PG yields specific information 
as to how the PG is distributed, and especially how likely it is that PG will be small (the RFI mim- 
ics the signal) or large (RFI is highly attenuated by the MF). Second, the designer can ask how the 
choice of a random ensemble for H affects the distribution of PG, and in particular if there is any 
choice that renders that distribution independent of r, achieving the goal of robust RFI rejection in 
a statistical sense. 

The IP suggests an immediate answer to the second question because it establishes that the energy 
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Figure 2: The cumulative distribution function Fpg(p) of processing gain PG for a CR signal H = X/||X||,a 
Gaussian CR X, and any r. It is plotted on a log-log plot (10 logjQ p vs logjQ F-p (p)) for K = 10" and n 6 [1, 8] . 
The distribution moves to the right by approximately 10 dB for every order-of-magnitude increase in K. 



of a Gaussian CR H is the same in any direction and thus the variance |H^r| does not depend on 
r. The only problem is that this ensemble does not incorporate the constraint that 1 1 H 1 1 = 1 . To 
address this, let X be a Gaussian CR random vector, and choose H = X/ 1 1 X 1 1 . For this choice of H, 
the distribution function of PG in ||9]l is ( Appendix B| 



Fpcip) = Ft {PG < p} = , p>l. (10) 

Notably FpQ{p) still does not depend on r and FpQ{p) decreases exponentialy with increasing K 
as illustrated in Figure |2] Any objective for the probability that PG < p for any p > can be met 
by choosing K sufficiently large. For example, if we require the probability that PG falls below 
40 dB to be less than 10"^°, then we must choose K > 230, 248 (e.g. T = 0.1 sec and B = 2.31 
MHz). 

Consider for example a jammer who deliberately tries to design a worst-case interferer and is 
aware that H was chosen as described previously, but is not aware of the specific outcome H. This 
jammer has no basis for choosing a specific r; that is, no choice is better or worse than any other. 
There is a duality: Choosing a CR signal H yields a PG distribution that is independent of r, just 
as a CR noise N results in a distribution (the same distribution!) that is independent of h. These 
two statements are compatible, as any randomly chosen H in the former is entirely suitable for the 
latter. For any other way of choosing H, the distribution of PG is dependent on r. That is, some 
interferers are rejected more than others, and robust immunity is lost. 

There are several conclusions. First, there are always values of r ~ H for which immunity is 
poor. In short, beware of impostors. However, the probability of this happening inadvertently 
can be made arbitrarily small by increasing K. Second, K can be chosen freely without penalty 
in noise sensitivity with a MF; large K is advantageous for RFI immunity and not deleterious to 
noise immunity. Third, robustness in RFI immunity requires that f be spread uniformly over its 
K DOF's, and this is achieved if H is Gaussian CR. Unlike noise immunity, the goal of robust RFI 
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immunity provides specific information about advantageous signal characteristics; namely, the 
signal should statistically resemble white Gaussian noiser] 




Choosing the signal from a CR ensemble in conjunction with the MF is an equilibrium strategy in 
the game theoretic sense. Any deviation from this strategy either degrades the detector sensitivity 
in noise, or it abandons robustness in RFI rejection, or both. Assuming the receiver designer fo- 
cuses on the limited information available to the transmitter designer, each will arrive at this same 
fundamental conclusion based on elementary probability theory. However, turning this principle 
into a concrete signal design introduces some ambiguity as discussed next. 

4. CONCRETE SIGNAL DESIGN 

Even acknowledging fundamental principles, obvious questions loom. How do we choose the 
signal parameters B and T? Given that the transmitter cannot communicate a chosen signal h to 
the receiver, can we make a credible argument that h can be guessed by the receiver, or that there 
are at least a limited number of options that can be searched? In addressing these issues, we draw 
upon extensive experience in designing terrestrial communication systems. 

4.1. Choosing T and B 

The RFI PG depends only on the time-bandwidth product K = BT, and thus offers no insight 
into how a chosen K should be partitioned between B and T. This design choice must rely on 
external considerations such as the transmitter's stellar scanning dwell time and desired data 
rate (which influence T), impairments in the interstellar medium (which influences B), and the 
receiver 's search strategy. Generally it can be asserted that T is a scarce resource but B is more am- 
ple. Another consideration is the RFI FG (Section |2.3| , which depends on the anticipated nature 
of the RFI. For persistent narrowband RFI (like other communication signals) it is more advanta- 
geous to increase B, but for bursty RFI (like impulsive noise from the electrical grid or motors) it is 
more advantageous to increase T. Finally, many current searches depend entirely on large FG by 
choosing K w 1 and thus perform poorly for RFI that mimics the signal's small time-bandwidth 
product but relatively well for other types of interferers. 

4.2. Pseudorandom signal 

In the random signal approach, there is no possibility of communicating h to the receiver, but the 
receiver may be able to guess an algorithm used to generate the signal. This suggests using a 
pseudorandom algorithm (PR) for generating a signal that is representative of the random ensemble. 
The criterion for "representative" is meeting standard statistical tests for a sequence that is CR. The 
significance or reliability of such tests increases with K. Although the PR algorithm mimics white 
Gaussian noise, the receiver searches for a unique signal rather than applying statistical tests to 



^ Since a Gaussian CR signal vector remains Gaussian CR after any transformation of basis, this conclusion is transparent 
to the basis. 
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the reception, and thus there is no extra source of confusion with noise when the signal is chosen 
randomly. 

4.3. QPSK signal 

To be successful in acquisition, the receiver must guess both the basis functions and the PR al- 
gorithm. Of course, in both cases the receiver can search multiple possibilities subject only to 
computational power (principally budget and technology) limitations and must suffer a penalty 
in PpA- (There is no possibility of searching over all possibilities, because the PpA penalty would be 
too large.) Based on our terrestrial experience, either the sampling theorem ||5) or Fourier ||6) basis 
seems a likely choice. When resolving the ambiguity of what PR generator to use, Occam's razor 
- look for the simplest solution that meets the requirements - is a good guiding principle [1]. The 
question then is whether the simplest possible PR generator, one that mimics a sequence of inde- 
pendent coin tosses, can be used. The answer is yes if an independent CR (rather than Gaussian 
CR) H can suffice. Fortunately, from Table |3] independent CR vectors obey the energy IP precisely 
and also approach Gaussian CR in distribution as K ^ oo. 

The approach used in terrestrial systems is to define a finite constellatioij^from which components 
of H are chosen randomly. The smallest complex-valued constellation for the components of H 
that can be independent CR and zero mean has four points (these properties are inconsistent with 
a two- or three-point constellation). This constellation must include the points (±1 ± f), where the 
real and imaginary components are independent and chosen by a fair coin toss. Suppose a binary 
PR generator statistically mimics a sequence of 2K independent coin tosses {c;, 1 < Z < 2K}, 
where c/ — ±1. Then the vector 

h=[ci+i- C2, C3 + f ■ C4, . . . , C2K-1 + i ■ C2k\ 

statistically mimics an independent CR random vector H, and it is the simplest choice. This is 
called a quaternary phase-shift key (QPSK) PR generator, and as expected the distribution of PG 
for this PR sequence and a Gaussian CR H are virtually indistinguishable as illustrated in Figure 

13 

4.4. The PR open-ended sequence property 

Another requirement on the PR algorithm is evident upon examining the basis expansions. The 
receiver's computational requirements can be dramatically reduced if a mismatch in DOF K be- 
tween transmitter and receiver is permitted, and this is permissible as long as the PR sequence 
{cjt, 1 < A: < oo} is open-ended and has an unambiguous initial state. If the bandwidth B remains 
fixed in ||5) and C]^ is mapped onto coordinates H with DOF K in the transmitter and M 7^ K in 
the receiver, then detection remains feasible with reduced sensitivity. If K > M then some trans- 
mitted energy is lost, if K < M some unnecessary noise is added into the signal level estimate, 
and maximum sensitivity occurs for M = K. The receiver must still search over different values 



This is similar to the constellation of PAM data symbol {Aj.}, but with a completely different motivation. 



14 




Figure 3: The distribution of RFI P in dB for an interferer consisting of a 64-QAM digital communication 
signal with random data symbols, one symbol per second, and a 50% excess bandwidth. The two cases 
shown both assume that T = 10 seconds, with B = 2 (J<C = 20) and B = 1000 {K = 10^). Shown are the results 
of one million trials of a PR QPSK signal using the RandomTrials[] function of Mathematica. The vertical 
axis is logjg of the fraction of trials for which P < p and the horizontal axis is p in dB. The solid curve is the 
theoretical prediction for a Gaussian CR signal. 

of B (or equivalently different sampling rates), but the value of T can be left open. This required 
search over B introduces an additional search parameter and increases the computational burden 
accordingly. Similar logic applies to Fourier basis except that the search is over T and it is B 
that is left open. 

Any discovery algorithm for an information-bearing signal must search over starting time and fc- 
In the discovery of an information-bearing signal (and spread spectrum is no exception) there is 
one additional dimension to search, which we call a dilation parameter. It can be chosen to be T 
(equivalent to the frequency sampling or symbol rate 1 / T) or B (equivalent to the time sampling 
rate 1 / B). The most sensitive detectors will search over both these dilation parameters. 

4.5. PR coin-toss generators 

Because of their explicit coordination, terrestrial systems do not require the open-ended sequence 
property and hence offer little direct guidance as to choice of a PR algorithm. Nevertheless it 
is straightforward to find PR algorithms with this property, including ones that rely on nothing 
more sophisticated than the geometry of the square or the circle. For example, common irrational 
numbers such as tt, \/2, and e expanded in base two satisfy the open-ended sequence property 
and have been observed to satisfy statistical tests for independent and identically distributed coin 
tosses [42]. If the receiver has to search over L different algorithms, an increase in £}, by log L is 
necessary to avoid an increase in overall PpA t23l- 

The number tt seems like a particularly interesting candidate, especially since it has been known 
and studied centuries and has been confirmed to pass rigorous statistical tests |43|. An example of 
one such test, the frequency of run lengths, is shown in Figure |4] A run of one +1 followed imme- 
diately by a -1 should occur close to one quarter of the time, two +l's followed by a -1 close to an 
eighth of the time, etc. Another test of randomness calculates the sample autocorrelation fimction. 
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Figure 4: Tests for independent coin-toss statistics for a two million bit binary expansion of n. On the left, the 
relative frequency of run lengths of -|-l's and — I's against the predicted probabilities. The disparities at small 
relative frequencies is explained by the statistically insignificant number of longer sequences. On the right, 
the magnitude of the autocorrelation sorted from largest to smallest is plotted with the unity autocorrelation 
at zero lag omitted. The predicted value of zero is satisfied within a few parts in a thousand. 



Since E[C,nC^] = ^ ■ 6m-n for any white CR vector, at large K the time average autocorrelation 
should approximate, 

Equivalently, the Fourier transform of cj^ should be approximately white (constant magnitude). 
By the IP, such signals have their energy spread relatively uniformly over both a time basis in 
f G [0, T], a frequency basis in / G [0, B], or indeed any basis. Such a signal appears noise-like and 
stationary, as illustrated by the n representation of Figures [5] and [6] 



5. CONCLUSIONS 



We have addressed the discovery of an information-bearing signal for interstellar communication 
in the presence of noise and RFl. Design for maximum detection sensitivity in noise offers no in- 
sight into the type of signal to transmit or seek. We have demonstrated that a signal in which each 
symbol-period pulse h{t) statistically matches a burst of bandlimited white noise is the only choice 
that achieves robust immunity to RFl (immunity not dependent on the specific in-band RFl wave- 
form) for matched filter detection. This immunity is additive to other forms of RFl mitigation, 
and works when RFl is overlapping in frequency, time, and space. A signal with this character- 
istic can be generated by a PR algorithm that we assert can be guessed by the receiver, such as a 
QPSK constellation in conjunction with the binary expansion of an irrational number like n. This 
observation informs both transmitter and receiver designers as to the desirable characteristics of a 
signal to transmit and to seek. 

A consequence of wider signal bandwidth that has not be addressed is its susceptibility to dis- 
persion and multipath distortion introduced in the ISM [|22| . While the same can be said of any 
information-bearing signal as the information rate is increased, spread spectrum increases the 
severity of this impairment. This is addressed elsewhere |23|, where it is shown that there is a 
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Figure 5: Plot of the real part of a PR signal with K = 100 based on the first 200 bits of the binary expansion 
of n (typically a much larger K would be used). By construction, the signal energy is spread uniformly over 
time, and is thus relatively immime to pulse-like RFI. 




Figure 6: Plot of the magnitude of the discrete Fourier transform of the signal in Figure |5| Reflecting the IP, 
the signal energy is spread relatively uniformly over frequency, and is thus relatively immime to narrowband 
RFI. The variation in energy with frequency is consistent with the expected complete randomness. 
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quantifiable no-dispersion region of {B, T, fc} for which these effects are insignificant. This ob- 
servation together with a desire to maximize K = BT suggests transmitting and searching for 
signals near the boundary of the no-dispersion region, dramatically narrowing the ambiguity in 
parameters {B, T, fc} for both transmitter and receiver 

Other practical issues are yet to be addressed. What are the consequences of relative motion be- 
tween transmitter, receiver, and ISM clouds? What about the conversion from continuous-time, 
and the related granularity in the search space over {fc,B} or {fc,T}7 What are the computa- 
tional requirements for such a search? How are these parameters related to a stellar scanning 
strategy? 
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Appendix A. INFORMATION-BEARING SIGNALS 

Assuming that the information to be conveyed is digital in nature Il30l , without loss of generality 
it can be represented as a sequence of bounded integers {4, — oo < k < 00} where k represents 
time and Ij. £ [1, M]. The radio electric field is generated by a process of modulation, turning {4} 
into a continuous-time waveform. In general form, and equivalently at baseband, this waveform 
is 

00 

y{t)= J2 Ki^-kT), (A.l) 

where the {hi (f ), 1 < Z < M} is a set of M waveforms distinctive in a way that allows the receiver 
to infer the corresponding 4- The raw information rate is K = (logj M) /T bits per second. 

A passband signal at carrier frequency fc can readily convey a complex-valued waveform y{t), 
which can effectively double or more than double the information rate by allowing a larger M. 
To see this, any real-valued passband signal x{t) with carrier frequency fc can be represented in 
terms of an equivalent complex-valued baseband y{t) as 

x{t) = V2 ■ U [y{t) ■ e'^^'f^'} . 

The \/2 factor equates the baseband and passband power or energy, and $R and S denote the real 
and imaginary parts of a complex number From 
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and y{t) can recovered from x{t) by lowpass filtering x{t) ■ e '^"M, 

Prominent special cases are orthogonal signaling, in which the {hni{t), 1 < m < M} are mutually 
orthogonal (uncorrelated), and the PAM of FSK and PSK Il28l are examples of these two signal 
classes. Although the choice of h{t) in PAM is emphasized in this paper, the results are readily 
extended to orthogonal waveforms. The results of this paper thus apply quite generally to any 
digital communication scheme conforming to I A.l\ (this includes virtually all digital transmission 
in use terrestrially), but not to analog modulation (which has been prominently applied to some 
METI experiments L8J). 



Appendix B. ISOTROPIC PRINCIPLE 



The statements in Table||will now be verified. A zero-mean complex-valued random vector in C 
can be written as V = A -|- f ■ B. A full set of second-order statistics include three K x K covariance 
matrices: E[AA^], E[AB"'"], and £[BB^] (x^ is the transpose of x). Equivalently we can use two 
complex-valued matrices (x^ is the conjugate transpose (Hermitian) of x) as [41 J 



R. 



and R7; 



VV' 



Rv is the conventional covariance matrix, and Rj, is a pseudo-covariance matrix. The equivalence can 
be verified by writing Rj, and Rj, in terms of A and B and solving for E [AA"*"], E [AB^], and E [BB""^] 
in terms of R^ and R^. In contrast to the real-valued case, the covariance matrix alone does not 
completely specify the second-order statistics. Further, it is readily shown that white CR in Table 



[2] is satisfied if and only if Ry = (7^1 and R^, = 0. 

Let U be a unitary K x K matrix (U''' U = U U''' = I) corresponding to a transformation from 
one orthonormal basis to another. For example U is the discrete Fourier transform (DFT) when 
converting from time to Fourier basis. If V is white CR, then so too is UV, as verified by 

(UV)(UV)^ =U(^E = (7^ -W^ = cr^ -l 

E [(UV)(UV)"'"j = E [uVV'^'u"'"] =\J (e [vv'^'j ) u""" = 0. 



Since UV is Gaussian when V is Gaussian, a Gaussian CR vector remains Gaussian CR after coor- 
dinate transformation. 

Consider the distribution of Zm and given by l|8|. By the IP, h''' Ni is Gaussian with variance 
cr-^ not dependent on h. For any K x K unitary U, e^'^ ■ U is also unitary and thus (since a unitary 
transformation is norm-preserving) 



/fft ■ Uh N2 



where by the IP since Ni is Gaussian CR then so too is N2 = e UNi. In particular choose any U 
such that Uh = [1 ... 0]^ (h is the first column of U'*', or equivalently the first basis vector). The 
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distribution of Zg does not depend on h because the distribution of N2 does not depend on U. For 
h = 

r K I 12 2A2l 

Pp. = Pr(ZH>A}=P,|E N_>^|, 

where the left side is and PpA can be calculated from the distribution function as shown 

in Figure [1] Similarly, Zm is PC^(2) and thus has the same distribution as Zg for X = 1. 

Our interest is in the random variable PG in ^ when h is a random vector H = X / 1 1 X 1 1 , and 
X is Gaussian CR. Let U be any unitary matrix whose first column is r. By the IP, V = X = 
[Vj, V2, . . . , Vk]^ is also Gaussian CR. Then from the norm preserving property PG is a function 
of r but its distribution is not, 

PQ_ W^f _ _ _ l|V||^ 

~ I X+ r I ^ ^ I X+ UU+ r I ^ ^ I V+ U+ r I ^ ^ \Vi\^ ' 

Letting W = Et=2 \ '^k\ > the distribution fimction is 

Pr{PG < p} = Pr|w < (p -1) ■ IVil^l = ^~Pr < (p - 1) ■ \ Yi\^^ f^N{w) dw 

where W is x^i'^^ ~ 2) and | Vi| ^ is ;t^(2). Evaluating the integral yields | |Tol l. 

The asymptotic property for an independent CR vector mentioned in Table |3] is significant for 
characterizing PG when H is independent CR but not Gaussian (e.g. a QPSK signal drawn from 
independent coin-tosses). The distribution of X^r will approach Gaussian as X 00 for almost all 
values of r since 

K-l 
k=0 

is a sum of independent random variables. Although the terms in the summation are not iden- 
tically distributed, the Lyapunov version of the Central Limit Theorem |44j asserts that 5?{X^r} 
and 3{X^r} approach Gaussian in distribution as K — )• 00 if at most a finite number of the 
{rj-, < A: < 00} are zero and the third moments of 5?{Xj-, < A: < 00} and S{Xjt, < A: < 00} 
obey certain upper bounds. For the QPSK PR generator of Section [43] these third moments are 
identically zero, so relatively rapid convergence to Gaussian is assured. Could any distribu- 
tions other than Gaussian result from the sum of independent random variables? Yes, the larger 
class of stable distributions |45|; however, the only stable distribution with finite variance (and 
£|X'''r| ^ = cr^ < 00) is the Gaussian. 

Appendix C. NOISE STATISTICS 

A noise vector N that originates as white Gaussian noise on the passband channel is Gaussian CR 
|[30 1 as long as /c > B. N is Gaussian because the conversion is linear. To demonstrate that N is also 
white CR, choose K orthonormal baseband basis functions {gk{t), ^ < k < K} on t & [0, T], and 
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assume real-valued white Gaussian noise M{t) at passband with autocorrelation E[M{t)M{s)] = 
No ■ S{t - s). Then 

E [NuNn,] = 2No e-'^^f^'gl{t)gl (t) dt = for /, > B (C.l) 

E [NkK] = 2No gl{t)gm{t) du = 2No ■ S^^,,, . 

Note that g"^ {t)gm (0 ^^11^ ™ frequency band / e [0, 2B], and i jCl^ is a Fourier transform evaluated 
at frequency 2/c, which evaluates to zero as long as 2fc > 26. 
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